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Resume : Nous presentons dans ce rapport la resolution, par une methode volumes finis, 
du systeme des equations de Saint- Venant avec termes sources topographiques sur des 
domaines ID. Avec une idee originale de Leroux Q], le systeme des equations est complete 
par une equation triviale sur la bathymetrie. Par un changement de variable, on elabore une 
formulation celerite-vitesse des equations que Ton linearise. Nous construisons ensuite un 
solveur de Riemann approche qui preserve la positivite de la celerite et qui assure la prise 
en compte des bancs-couvrants-decouvrants. Enfin, des applications numeriques sur des cas 
tests sont presentees. 
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A VFRoe scheme for ID shallow water flows : 
wetting and drying simulation 

Abstract: A finite- volume method for the one-dimensional shallow- water equations includ- 
ing topographic source terms is presented. Exploiting an original idea by Leroux j^, the 
system of partial-differential equations is completed by a trivial equation for the bathymetry. 
By applying a change of variable, the system is given a celerity-speed formulation, and lin- 
earized. As a result, an approximate Riemann solver preserving the positivity of the celerity 
can be constructed, permitting wetting and drying flow simulations to be performed. Finally, 
the simulation of numerical test cases is presented. 

Key-words: Shallow water equations, finite volumes, Riemann solver, positivity preserving 
scheme, wetting and drying fiows 
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1 Introduction 

La ville de Cotonou est dans une cuvette ceinturee de plans d'eau et traversee du nord 
au sud par un chenal (la Lagune de Cotonou) de 4000 m de long sur 350 m de large, reliant 
le lac Nokoue a I'ocean Atlantique (FigQ. 

Les quelques ouvrages d'assainissement dent dispose la ville sont generalement en bordure 
des routes et, pour la plupart, connectes au chenal. Malheureusement, ce systeme de drainage 
au lieu d'aider a I'evacuation des eaux pluviales, sert plutot de vecteur a I'invasion de la 
ville par les eaux de crue du chenal. 



Les ecoulements d'un tel reseau d'adduction d'eau sont des ecoulements a surface libre, 
en eaux peu profondes (shallow water) et sont alors regis par un systeme bidimensionnel des 
equations de Saint- Venant. Ces equations sont obtenues a partir des equations de Navier- 
Stokes pour un fluide incompressible en faisant I'hypothese de pression hydrostatique, de 
vitesses uniformes suivant la verticale, d'un fond et d'une surface libre impermeables. On les 
emploie dans des domaines aussi divers que la protection de I'environnement, le calcul des 
marees et des ondes de tempete, la sedimentologie, la simulation des ondes de submersion, 
I'etude des crues, etc. 

Dans ce rapport, comme premiere etape dans la realisation d'une simulation numerique 
realiste, nous proposons une etude du systeme unidimensionnel des equations de Saint- 
Venant avec topographic. Suivant I'idee developpee dans plEilij, nous utilisons une meth- 
ode de type Volumes-Finis-Roe (VFRoe) se basant sur une formulation celerite-vitesse des 




Fig. 1 Chenal de Cotonou 
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equations. Une telle formulation pent conduire a une celerite negative de I'etat intermediaire 
du probleme de Riemann y decoulant. Notre contribution est la construction d'un solveur 
de Riemann approche avec un choix des vitesses d'ondes garantissant la preservation de la 
positivite de la celerite de I'etat intermediaire. 

Ensuite, grace aux resultats de jSHHI nous adoptons le solveur de Riemann approche afin 
de prendre en compte le cas des bancs-couvrants-decouvrants (zones noyees /non-noyees) . Le 
traitement de ce type de cas est essentiel en vue de la simulation numerique des problemes 
d'intrusion d'eau dans la ville de Cotonou. 

Enfin, pour tester la robustesse de notre schema, dans un premier temps nous reprenons 
des cas-tests effectues dans ^ puis dans un second temps deux derniers cas-tests sont presen- 
tes pour tester le debordement de I'eau d'un canal (bancs-couvrants-decouvrants). 



2 Le modele mathematique 



En absence des forces de frottement, le systeme ID des equations de Saint-Venant est 
donne par : 



( dh d{hu) 
dt dx 

dihu) 9 /, 2 1 ,2 



= 



, da 



,{x,t) e 



(1) 



u{x,t) est la Vitesse de I'eau, h{x,t) la hauteur d'eau, a{x) la hauteur de la topographic du 

sol, h + a I'elevation de la surface libre de I'eau (a et h + a sont prises par rapport a un plan 

de reference), g I'acceleration gravitationnelle. 

La bathymetrie a etant independante du temps, on complete le systeme (QJ par I'equation 

da ^ ,,_ , . 

— = 1 et on obtient : 
ot 



dh 


dq 


'dt " 


dx 


dq 


d 


-dt' 


dx 


da 




-dt 





1 



gh^ + gh 



da 
dx 



,(a;,t)eMx 



(2) 



ou q — hu represente le debit d'eau. 



RR n° 0123456789 



6 



A.W.BELLO, A.GOUDJO, C.GOUDJO, H.GUILLARD & J.-A.DESIDERI 



3 Integration numerique 

Nous presentons dans cette section I'integration des equations par utilisation des tech- 
niques de volumes finis et de solveurs de Riemann. 

3.1 Le maillage 

On considere un ecoulement unidimensionnel sur une longueur L. 

L'intervale [0,L] est subdivise en N segments de meme amplitude Ax — —. On obtient 
alors une suite de points (a;j")jej={o....,JV} definis par : Xj = jAx. On pose 



Co 



xo;xi 



j€ {1,...,7V-1} 



(3) 



^ 



Fig. 2 Discretisation spatiale 

Pour la discretisation temporelle, on se donne un pas de temps Ai et une suite d'instants 
t" = nAt, n > 0. 

3.2 Un schema de type Godunov 

La topographic a du probleme Q est approchee par des valeurs discretes Oj, j & J 



Ax 



a{x)dx. 



(4) 



En posant : 
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aA{x) ^^^ajXjix), (5) 
avec Xj{x) = I si X £ Cj et Xj{x) = sinon, le systeme Q est approche par : 



r dh 



dt dx 

dq ^ d f q^ I , 2\ . 1 '^^^ 



dt dx \ h ^ 2^ ) ^ ^ dx 



daA 



,(a;,i)eKx 



(6) 



I dt 



Posons w = q .La solution w du probleme (O est approchee par des valeurs 
discretes w", j G J, n g N : 



^7 ^ l^^i-^ndx. 

duA 



(7) 



Puisque, sur chaque cellule Cj, — — = 0, I'integration de © sur Cj x nous 



donne 



dt dx 



{w{x,t"+')-w{x,t")}dx + J^^ [Fiw{xJ^,,t))-F{w{x+_,,t))] 



dt = 



, 1 \ 

on F{'w) =1 ^ + hgh^ \ ] x et a;+ designant respectivement les limites a gauche et a 

, J 
droite de x. 

La premiere integrale de cette derniere relation est approchee en utilisant H restera 
alors I'approximation numerique de la deuxieme integrale. 

La methode des volumes finis repose sur le fait qu'a tout instant, la solution w est con- 
stante par cellule. Ainsi, partant de la solution t") a I'instant t", le calcul de w{xj_^i, t) 
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et w{x'^_i, t) pour t e est donne par la resolution du probleme de Riemann suiv- 



ant : 



dh 
dt 



dq 

dx 



dq d f q^ 1,2 
-d-t^d-x[-h^2'^ 



da/ 



gh 



daA 
dx 



dt 



= 



= 



= 



(8) 



f Wl, si a; < x,, i 
w(x,t'^) = \ . 

1^ wji, SI a; > Xjj^i 

ou Wl — wJ et Wfi — w".^^. 

Designons par wJ^-^^2i^/^i'^L,wji) la solution autosimilaire de ©■ 
On definit deux flux 

F:^;-,^F{w''I^,^,{Q-;wl,wr)) et F^^+, = F{w^^,/^{0+;wl,wr)) 

correspondant a chaque cote de I'interface de Xj^i. On obtient alors le schema numerique 
suivant : 



w 



n+l 



Ax \ 3+1 



' 1 + 



(9) 



4 Resolution du probleme de Riemann 

II vient de ce qui precede que I'implementation numerique des equations de Saint Venant 
reside en la resolution, a chaque interface, du probleme de Riemann l|Hl ; puisque c'est cette 
resolution qui fournira les flux numeriques a utiliser. 

En posant c = 
equivalent a : 



gh, a chaque interface a; = et en dehors des zones seches, Q est 
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( d{2c) d{2c) du 
dt dx dx 

du ^ d{2c) ^ du ^ 
dt ^ dx ^ dx 



daA 
' dx 



daA 
"df 



w{x,0) = 



wl, si x < 
wr, si a; > 



(10) 



En posant Y{w) 




on obtient finalement : 



dY ,,,,,dY „ 



Y{x,0) 



si a; < 
si a; > 



avec A{Y) 



u c 
c u g 




(11) 



4.1 Linearisation du probleme de Riemann 

A chaque interface, on resout le systeme lineaire suivant 



dY ,,,y,dY „ 



Y{x,0)=Yo{x) 




(12) 



ou y = Y{Yl, Yji) est un etat moyen dependant uniquement des etats Yl et Yh ; F doit en 
plus verifier la condition de consistance Y{Y*,Y*) = Y*. Un choix possible est done : 



Y = 



Yr + Yl 



(13) 
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4.2 Solveur de Riemann 

A{Y) admet trois valeurs propres qui sont : 

Ao = 0, Ai = u — c, A2 = w + c, (14) 

Sous I'hypothese que |m| 7^ c ct cn dehors dos zones seches, les trois valeurs propres sont 
deux a deux distinctes et A{Y) est diagonisable. 

A chaque valeur propre Afc, nous associons un vecteur propre a droite G : 

A{Y)ru = XkTk (15) 

et un vecteur propre a gauche : 

%A{Y) = Xk% (16) 

Proposition 4.1 

%rk = 0, Vj ^ k. 



Preuve 

t 



IjA = XjHj =^ {%A) rk = (a/Zj) rfc 

=> Hj {Ark) = XjHjVk 
^IjiXkVk) = XjHjVk 
=^ (Aft - Aj) HjVk = 



t 



IjVk = 0, Vi ^ k. 



Remarque 4.1 



''3' K — ''J 

• Sir est un vecteur propre de la valeur propre X, alors pour tout a ^ 0, ar est egalement 
un vecteur propre de X. On peut done choisir les bases (ro,ri,r2) et (lo,li,l2) telles que 
= 1, Vj = 0,l,2. 

On peut prendre : 
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et 



• ^0 



/ 3(A2-Ai) \ 
2 * 



g(A2+Ai) 

2 

V A2A1 J 



,ri=\ 1 I , r2 = I 1 







/ \ 



V A2A1 / 



V 2A, / 







/ '2 \ 



V 2A2 / 



lfT2ll etant strictement hyperbolique, alors (ro,ri,r2) est une base de R^. 
En designant par Y* la solution de ltT2|l . on pent alors ecrire : 
Y*{x, t) — aQ{x, t)ro + ai{x, t)ri + a2{x, t)r2 et 
Y^{x) = a[5(x)ro + a\{x)ri + a2{x)r2 



dY » dY 



i=0 
1 



i=0 
2 



i=0 



D'apres la DroDositior l4.lL on a : 

1 



vc\{x) 



-Yo*{x)-k^Yo*ix)-h 



af = Yl- k, si a: < 



af^ — Yff ■ k, si a; > 



Mais pour chaque i, I'equation d'advection 

dt ' dx 

0/ \ / (^i^ si X < 
* ^ ' a", SI a; > 



(17) 
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admet pour unique solution : ai{x,t) = 



af, si f < A, 



af , si f > Ai 



Ainsi la solution du probleme de Riemann Ijl2|l est 

x/t>\k 
2:/t<Afc 



4.3 Solutions du probleme de Riemann 

Nous presentons dans ce paragraphe les trois cas possibles de solution du probleme de 
Riemann II12II dans le plan {x,t). On designe Y*{0~,t) et F*(0+,t) par Yj* et Y* respec- 
tivement. 



u + c i 


^ t 




^ ^ \ 
u — c \ ' 


1 J. 




Yl X\ 




Yr 


(a) 


^ 







I u — c 






Y* / 


Yl 




Yr 

Myl ^ 



(b) 



u + c 



Fig. 3 Ecoulement torrentiel 



(a) 



y; = y; - [a^^ - 



Y: = Yr 



;(b) 



Y* = Yr 



Y* = YC + (a« - «^)ro 




u + c 



y; = y^? - (af - a^)r2 



Fig. 4 Ecoulement fluvial 
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4.4 Correction entropique 

Le solveur de Roe peut conduire a une celerite negative faisant ainsi exploser Timple- 
mentation numerique. 

3 5 

Par exemple, prenons un ecoulement sur fond plat (a(x) = 0) avec ul = ~2^' ~ 2^' 

1 - 3 

On a alors et djl) : Ai = --c < ; A2 = -c > 0. 

L'ecoulement local est done fluvial. Determinons la composante 2c;* de Y* (voir FIGlJl : 
2c;* — 2cl + (af — ai )r\ (avec r\, premiere composante de ri) 



2cL 



-^(2cfl - 2cl) + ^(ui?^ - Ul) 



2CL 



1,0 o ^ , 1 A 2cfl + 2cL 



= -{cr + cl)<0 

4.4.1 Celerite de I'etat intermediaire 

En reprenant les trois conflgurations de la section- IOl on obtient 
Cas (a) : 



Cas (b) 



Cas (c) 



2c*i = 2cR 
2cl = 2cR 



g{aR - aL){M - Ai) 
2A2A1 




g{aR - aL){M - Ai) 
2A2A1 



(18) 



(19) 



2cr 



- ( 2{cR + cl) - {uR - Ul) 



2c* ^ -[ 2{cr + Cl) - {ur - Ul) - 



gja-R ~ o-l) 
Ai 

gjaR - cil) 
A, 



(20) 



RR n° 0123456789 



14 



A.W.BELLO, A.GOUDJO, C.GOUDJO, H.GUILLARD & J.-A.DESIDERI 



4.4.2 Ecoulement sur fond plat 



Proposition 4.2 Si 



ur-ul < 2{cr + cl) 



(21) 



alors I'etat intermediaire Y* admet une celerite positive. 
Preuve 

La preuve est immediate; en effet, ici ~ — 0. Par suite, la celerite de I'etat 
intermediaire Y* est positive d 'anres [181191 et [2(11 



Si la condition lf2T|l est violee, on pose [H : c* = = 0. 

4.4.3 Ecoulement avec topographie 

On choisit alors d'agir sur les vitesses d'ondes A de I'etat intermediaire Y* ; on notera A 
les nouvelles vitesses de Y* qui seront utilisees dans I'implementation numerique. 

Proposition 4.3 On suppose ur — ul < 2{cr + cl) ; si Xi < < A2, le choix 



□ 




et 




permet Vohtention d'un etat intermediaire Y* a celerite positive. 



Preuve 



La preuve est immediate d 'anres [2(11 



□ 



Si la condition lf2T|l est violee, on pose : 

' 4Ai , si flfl > GL et < 

c*=0 \ 



cf =0 



c* = - 



gjaR — ol) , si afl; < 
4A2 
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Proposition 4.4 Si X2 < 0, le choix 
' Ai = Ai 

A2 = min A2, 



g(aR-aL)\\ \ , si ur > ql | _ ' ^ "'^ 

4c_rAi - g{aR - ql) I 



permet I'obtention d'un etat intermediaire Y* a celerite positive. 
Preuve 

La preuve est immediate d'a.nreslTsl 
Proposition 4.5 Si Ai > 0, le choix 



r \ \ I \ (\ g(afl - aL)A2 ^ 

Ai = Ai . J Ai = max Ai, 

1 A ^ , SI QR > aL et < y 4clA2 - g{aR - ul) J > si ur < cl 



A2 — A2 

permet I'obtention d'un etat intermediaire Y* a celerite positive 



Preuve 

La preuve est immediate d'anreslTol 

□ 



4.5 Traitement des bancs-couvrants-decouvrants 

On se situe ici dans le cas ou la hauteur de I'etat wr ou celle de wl est nulle. Pour la 
determination de I'etat intermediaire Y* , nous utilisons les resultats de bases sur une 
nouvelle estimation des vitesses d'onde Ai et A2 : Une onde de detente est generee du cote 
on la hauteur d'eau est nulle jHEI- 
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= ul-cl 
: ul + 2cl 




= ul- 2cl 
Ul + cl 



Fig. 5 - bancs-couvrants-decouvrants 



5 Resultats numeriques 

5.1 Ecoulements noyes 

La longueur du domaine dc calcul est L = 25m, le nombre de points du maillage N = 
1000 ; la duree totale d'observation est T = 1.2s ; les conditions initiales sont : 





< X < 12.5 


12. .5 < X < 25 




3 


4 


u{x, t = 0) 
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5.1.1 Testl : Ecoulement sur fond plat 

a{x) = 




Fig. 6 - Testl (t=Os) 




Fig. 7 - Testl (t=1.2s) 
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5.1.3 Tests : Ecoulement avec topographie 

a{x) = 




Fig. 10 - Tests (t=Os) 



Eau 



Fig. 11 - Tests (t=1.2s) 
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5.2 Bancs-couvrants-decouvrants 

La longueur du domaine de calcul est L = 25m, le nombre de points de maillage N = 500 ; 
la duree totals d'observation est T ; les conditions initiales sont : 

Test 4 





Q<x< 12.5 


12.5 < a; < 25 


h{x,t^ 0) 


1.5 





u{x, t = Q) 









Test 5 





< X < 7.5 


7.5 < .X < 10.5 


10.5 < a; < 14.5 


14.5 < a; < 25 


h{x,t = 0) 


1 


0.1 -a{x) 





0.1 - a{x) 


u{x, t = 0) 
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5.2.2 Tests : Bancs-couvrants-decouvrants 

a{x) (' 




IS 



•1 -- 



Fig. 14 - Test5 (t=Os) 









Eau 




n 






Eau 



IS JJ 



Fig. 15 - Test5 (t=3s) 
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6 Conclusion 

Exploitant les resultats de [UEl^, nous avons presente une formulation celerite-vitesse 
des equations de Saint-Venant. La methode VFRoe conduit a la resolution, par interface du 
maillage spatial, d'un probleme de Riemann. Notre contribution a ete, en agissant sur les 
vitesses d'ondes, la construction d'un solveur de Riemann approche garantissant la preser- 
vation de la positivite de la celerite de I'etat intermediaire dans la resolution du probleme 
de Riemann. L 'exploitation des resultats de nous a ensuite permis d'adapter le solveur 
afin de prendre en compte la simulation numerique des bancs-couvrants-decouvrants. 

Le solveur de Riemann approche presente dans ce rapport permet d'effectuer des simula- 
tions numeriques des equations de Saint-Venant, meme en presence de hauteurs d'eau nulles 
(recouvrement, decouvrement). La prochaine etape de notre travail consistera a I'extension 
a des tests bidimensionnels. 



RR n° 0123456789 



24 



A.W.BELLO, A.GOUDJO, C.GOUDJO, H.GUILLARD & J.-A.DESIDERI 



References 

[1] Ashwin CHINNAYYA, Alain- Yves LEROUX, and Nicolas SEGUIN. A well-balanced 
numerical scheme for the approximation of shallow-water equations with topography : 
the resonnance phenomenon. International Journal on Finite Volume , 1-1 :l-33, 2004. 

[2] Thierry GALLOUET, Jean-Marc HERARD, and Nicolas SEGUIN. Some approximate 
godunov schemes to compute shallow-water equations with topography. Computer and 
Fluid , 2001. 

[3] Thierry GALLOUET, Jean-Marc HERARD, and Nicolas SEGUIN. Some recent finite 
volume schemes to compute euler equations using real gas eos. International Journal 
For Numerical Methods In Fluids , 39 :1073-1138, 2002. 

[4] Thierry BUFFARD, Thierry GALLOUET, and Jean-Marc HERARD. Un schema sim- 
ple pour les equations de Saint-venant. http ://www. cmi.univ-mrs.fr/~gallouet/Publi/tg-art981- 
svt.pdf. 

[5] David L. GEORGE. Numerical Approximation of the Nonlinear Shallow Water 
Equations with Topography and Dry Beds : A Godunov- Type Scheme . PhD thesis. 
University of Washington, 2004. 

[6] Randall J. LEVEQUE. Finite volume methods for hyperbolic problems , pages 203-252. 
Cambridge University Press, 2002. 

[7] Edwige GODLEWSKI and Pierre- Arnaud RAVIARD. Numerical Approximation of 
Hyperbolic Systemcs of Conservation Laws , pages 1-98. Springer, New- York, 1996. 

[8] Herve GUILLARD and Remi ABGRALL. Modelisation numerique des fluides 
compressibles . Gauthiers-villars, 2001. 

[9] A. BERMUDEZ, A. DERVIEUX, J.-A. DESIDERI, and M.E. VAZQUEZ CENDON. 
Upwind schemes for the two-dimentional shallow water equations with variable depth 
using unstructured meshes. Comput. Methods Appli. Mech. Eng. , 1998. 

[10] Randall J. LEVEQUE and David L. GEORGE. Hight-resolution finite vol- 
ume methods for the shallow water equations with bathymetry and dry states, 
http ://www. amath.washington.edu/~rjl/pubs/catalina04/catalina.pdf, February 2006. 

[11] Olivier THUAL. Propagation des ondes de surface, articles pedagogiques, 2004. 

[12] Philippe HELLUY. Approximation des systemes de lois de conservations hyperboliques. 

http ://helluy.univ-tln.fr/PHYP/phyp.pdf. 



INRIA 




Unite de recherche INRIA Sophia Antipolis 
2004, route des Lucioles - BP 93 - 06902 Sophia Antipolis Cedex (France) 

Unite de recherche INRIA Futurs : Pare Club Orsay Universite - ZAC des Vignes 
4, rue Jacques Monod - 91893 ORSAY Cedex (France) 
Unite de recherche INRIA Lorraine : LORIA, Technopole de Nancy-Brabois - Campus scientifique 
615, rue du Jardin Botanique - BP 101 - 54602 Villers-les-Nancy Cedex (France) 
Unite de recherche INRIA Rennes : IRISA, Campus universitaire de Beaulieu - 35042 Reimes Cedex (France) 
Unite de recherche INRIA Rhone-Alpes : 655, avenue de I'Europe - 38334 Montbonnot Saint-Ismier (France) 
Unite de recherche INRIA Rocquencourt : Domaine de Voluceau - Rocquencourt - BP 105 -78153 Le Chesnay Cedex (France) 



Editeur 

INRIA - Domaine de Voluceau - Rocquencourt, BP 105 - 78153 Le Chesnay Cedex (France) 

http: / /www. inria.fr 
ISSN 0249-6399 



I A 

F UTURS 



FLIA 

LORRAINE 



INKIA 

RENNES 



INKIA 

RHONE-ALPES 



INKIA 

ROCQUENCOURT 




INFIIA 



WINRIA 



apport 

technique 



